Field-Free Synthetic-Ferromagnet Spin Torque Oscillator 
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We study the magnetization dynamics of spin valve structures with a free composite synthetic 
ferromagnet (SyF) that consists of two ferromagnetic layers coupled through a normal metal spacer. 
A ferromagnetically coupled SyF can be excited into dynamical precessional states by an applied 
current without external magnetic fields. We analytically determine the stability of these states 
in the space spanned by the current density and SyF interlayer exchange coupling. Numerical 
simulations confirm our analytical results. 
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The transfer of angular momentum between the mag- 
netic layers of current-driven spin valves (spin-transfer 
torque) has not been so long ago predicted^i^ and ex- 
perimentally confirmedi^i^ The imphed efhcient electrical 
control of magnetizations motivated the pursue of new 
research directions. When the current density exceeds 
a critical value, the spin-transfer torque can switch the 
magnetization to a different static configuration without 
the necessity of applied magnetic fields, which makes it 
attractive for next generation Magnetoresistive Random 
Access Memory (MRAM) application^^j^^^ Under an ex- 
ternal magnetic field, the spin-transfer torque can also 
drive the magnetization into sustainable coherent oscil- 
lations spanning a wide frequency range from a few MHz 
to several hundred GHzi^i^^— High frequency magnetic 
oscillations generate a coherent microwave voltage signal 
through the Giant Magnetoresistance (GMR) in metallic 
spin valves or through the Tunneling Magnetoresistance 
(TMR) in magnetic tunnel junctions (MTJs). This effect 
can be used in so-called spin-torque oscillators (STO), 
which has many advantages including wide tunability,!^ 
very high modulation rates ;^'^'^^ compact device size, and 
high compatibility with standard CMOS processesJ^ii^ 
Thus STO is appealing for high frequency microwave ap- 
plications including microwave emitters, modulators and 
detectorsiii However, the necessity of an applied mag- 
netic field up to ^1 Tesla has greatly limited the poten- 
tial of these STOs for microwave generation and wireless 
communication applications. Recently, various solutions 
have been proposed to enable zero-field operation, viz. 
STO with a perpendicularly magnetized fixed^^ or spin 
valves with out-of-plane magnetized free layer^^^i^ mag- 
netic vortex oscillators)^"— wavy-torque STO by judi- 
cially choosing free and fixed layer materials with differ- 
ent spin diffusion lengths?^ and a tilted magnetization 
of the fixed layer with respect to the film plane.—"— 

Recently, synthetic ferromagnets (SyFs) composed of 
two ferromagnetic layers separated by a very thin non- 
magnetic spacer have been used to replace the free layer 



of a spin valve or MTJ<^"— SyF based spintronic devices 
have the advantage of higher thermal stability, smaller 
stray magnetic fields, faster switching speed and reduced 
threshold switching current as compared to single fer- 
romagnetic free layers;^"— Klein et alM- predicted that 
an anti-ferromagnetically coupled SyF layer with uncom- 
pensated magnetization can generate microwave oscilla- 
tions at zero applied magnetic field. 

Here we predict that a ferromagnetically coupled SyF 
can also be driven into dynamical precessional states, 
which, however, are surrounded in parameter space by 
static canted states with non-coUinear magnetizations. 
We use an analytical approach to determine the stability 
regimes of the SyF system and confirm results by numer- 
ical simulations. 

We study a spin torque nanodevice with synthetic fer- 
romagnetic free layers as shown in Fig. [T] The left 
ferromagnetic film forms the fixed polarizer with mag- 
netization mo||z, and the SyF consists of two ferro- 
magnetic layers FMi and FM2 of thickness di,2 with a 
paramagnetic spacer. The unit vectors describing the 
magnetization orientation are mi for FMi and 10.2 for 
FM2. For simplicity, we assume that the SyF layers 
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FIG. 1. A spin valve structure with an SyF free layer, where 
FMo is the fixed layer and FMi,2 layers are (anti-) ferromag- 
netically coupled. 



2 



are made of the same materials with identical satura- 
tion magnetization Mg. The exchange coupling strength 
reads Ec — —JSmi m.2, where J and S are the cou- 
pling energy per unit area and the cross section area of 
the sample, respectively. This corresponds to an effec- 
tive coupling field — Jmi/{fioMsdi), where i = 1,2 
and i = 3 — i, fiQ is the vacuum magnetic susceptibil- 
ity, mi and ni2 can be parallel or anti-parallel at zero 
applied field, corresponding to the non-local Ruderman- 
Kittel-Kasuya-Yoshida (RKKY) exchange ferromagnetic 
(J > 0) or antiferromagnetic (J < 0) coupling, respec- 
tively. The spacer between FMq and FMi is presumed 
thick enough that the RKKY coupling with the fixed 
layer is negligibly small. Although the dynamic dipo- 
lar coupling may be responsible for the apparent reduc- 
tion of static magnetization^ or linewidth of the current- 
induced spin wave mode^, it is estimated to be much 
smaller for our case compared to the shape anisotropy 
field and the other fields due to current-induced spin 
torque and interlayer exchange coupling and therefore 
disregarded^^. 

Let Pq^i be the spin current polarization by mo,i such 
that the spin current density in the two spacers are PqJ 
and Pij with j the electric current density. The corre- 
sponding spin-transfer torques on mi and m.2 are given 
by the projections: 
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with 7 the gyromagnetic ratio and Pq (Pi) are in general 
functions of the angle 9 = Z(mo,mi) (Z(mi, m2)).— 

Spin pumping causes enhanced damping in a ferromag- 
netic layer by emitting spin current into the adjacent 
non-magnetic layers4^ This emitted spin pumping cur- 
rent can exert a torque on the second layer. Disregarding 
the backflow and diffusion in the spacer layer, the torque 
density acting on m^ due to spin pumping from mq can 
be written as 



Nspi I3mjxmj - [(/Sm^ xmj)-mi]mj 



(2) 



where /3 is the effective enhanced damping due to spin 
pumping. It has been shown that Eq. ([2]) gives rise to a 
dynamic exchange interaction that can induce synchro- 
nization of the magnetization dynamics in two neighbor- 
ing ferromagnetic layers even for wide spacers.-^ In the 
results below we fully include the spin pumping. How- 
ever, in contrast to multilayers excited by microwaves,^'* 
we observe here only small corrections demonstrating the 
dominance of charge current-induced torques. 

The dynamics is described by the coupled Landau- 
Lifshitz-Gilbert-Slonczewski (LLGS) equationsj^i^ 



where a is the sum of the intrinsic Gilbert and the 
spin pumping induced dampingJi The effective mag- 
netic fields Hi consist of shape anisotropy and RKKY 
exchange coupling and can be written as. 
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For simplicity, we consider di = d2 = d (equal magne- 
tization) for the rest of the paper unless otherwise speci- 
fied. We linearize Eq. ([3]) in the vicinity of four collinear 
equilibrium states, i.e. ttj ti; 4-tj ii; and assume m.^ — 
XiZ + Ui with Xi = ± and Ui denoting the small trans- 
verse magnetization component. After the linearization 
and the Fourier transform Ui{t) = Jui(w)e~"^*d(u;/27r, 
Eq. ([3]) becomes 
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with ujQ = 2-fKu/^J.aMs, usj = jJ/fioMsd, ujj = 
{h/2e){jj/iJ,QMsd). The frequency of the normal modes 
are given by the eigenvalues oi W — —A^^V : fli and 
ri2- When any of the Im r2i_2 > 0, the system is unsta- 
ble, implying that an infinitesimal perturbation will lead 
to magnetization dynamics with amplitudes that initially 
increase exponentially in time. 

The above results allow us to calculate the stability re- 
gions for the tt) til itj a phases in the space of typical 
experimental parameters: angle- independent Pi = P2 = 
P = 0.5, d ^ 3 nm, Ku = 8x10"* J/m^, j^lO^ A/cm^ 
and J^l mj/m^i^i^ To analytically construct the sta- 
bility diagram as shown in the top- left panel of Fig.^a), 
we first calculate the eigenvalues for each given set of 
[j,J] as given by Eq. ([5]). Then we determine whether 
any of the four collinear static states (different combi- 
nations of [Ai,A2]) is stable or not. For example, both 
the imaginary part of the eigenvalues of ft configuration 
[Ai — -1-1, A2 = -1-1] are negative when j<0. Therefore 
tt is stable in the blue region. In this way, we quickly 
map the parameter space for any given set of [j, J] and 
construct the entire stability diagram consisting of four 
collinear magnetization configurations. The spin torque 
drives the SyF to the parallel tt configuration for nega- 
tive currents j. For positive currents, the 4,t configura- 
tion is preferred. These results can be understood from 
Eq. ([T]). In a small region the antiparallel fi state exists 
for negative J and small j {i.e. in the vicinity of the neg- 
ative vertical axis but not visible in the figure due to the 
scale). Although it seems that the H state also occupies 





5 10 20 40 5 10 




5 10 20 40 20 40 

t (ns) t (ns) t (ns) 



FIG. 2. (Color online) Dynamical phase diagram in the pa- 
rameter space of currents and RKKY coupling strengths, (a) 
Phase diagram calculated analytically by Eq. ([5]) ; none of the 
four states tt; tii -Ul is stable in the white region, (b) 
Phase diagram calculated by numerically solving the LLGS 
Eq. ([3)|. The purple are the STO phase, and the white one 
the canted state, (c) The time evolution of the polar angles 
6\^2 ~ Z(mi,2,z) at the six different points indicated in the 
phase diagram. In the third subfigure of (c), the solid and 
dashed lines correspond to different sets of initial conditions. 



the fourth quadrant (j' > 0, J < 0), this triangular region 
is hysteretic, i.e. 4-i and |t m^Y both appear depending 
on the history. 

Most importantly, there is a white/purple region in 
which none of the four static collinear states is stable, 
therefore it must be either in a dynamical STO or static 
canted state. To leading order of a , we find from Eq. ^ 
an approximate boundary for the white region: 

upper: ujj = ujj, (7a) 

lower: w/ = iAwq + oj^ - 2ujo + a—, (7b) 
V cjj 

which is plotted as the black dashed lines in Fig. [2ra,b), 
matching the numerically obtained boundaries almost ex- 
actly. Eq. ([7]) is calculated from the eigenvalue anal- 
ysis based on Eq. ([5]) with perturbation from the four 
static collinear states. This method is equivalent to that 
used by Bazaliy et aZ."*^ A fully analytical solution for the 
boundary between STO and static canted phase turned 
out to be intractable, due to the complexity of Eq. ([5]) 
for non-collinear states. 

We now present numerical solutions of the LLGS 
Eq. ([3]) including damping, spin torque and RKKY cou- 
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FIG. 3. (Color online) Power spectrum for mo-mi (left) and 
ini-m2 (right) as a function of current density j and frequency 
/ at J = 0.25 mJ/m^, corresponding to the black line in the 
top right panel of Fig. [21 



pling. We summarize the dynamics of the coupled mi 
and m.2 in Fig. [Hb), in which we confirm the phase 
boundaries in the analytical analysis in Fig. HJa). In 
addition, we can now map the STO phase by the pur- 
ple color. The rest of the white region consists of static 
canted states. In Fig. ^c) we show the six different SyF 
configurations that may exist depending on the current 
and RKKY coupling strength. Point 5 corresponds to an 
STO state, in which both mi^2 are undergoing large an- 
gle precessions, which result in a large magnetoresistance 
oscillations attractive for applications. 

For the STO phase, we study the power spectrum 
of the magnetoresistance due to the magnetization os- 
cillation of mi^2, which is approximated by R{t) = 
i?o + A_Rimo-mi + A_R2mi-ni2. Fig. [3] shows the Fourier 
transform of mo-mi (left) and mi-m2 (right) as a func- 
tion of current density j at J = 0.25 mJ/m^ (correspond- 
ing to the black line in Fig.[21[b)). The clear higher order 
harmonic modes are evidence of the non-linearities in the 
STO dynamics. Fig. [3] also demonstrate that the oscil- 
lation frequency of the device can be continuously tuned 
by the current at zero applied magnetic field and thus 
potentially be utilized for nano-scale microwave applica- 
tions. It should be noted that the frequency range can be 
further tuned by tens of GHz by adopting a larger or 
taking into account the easy-plane anisotropy field (de- 
magnetization field) . 

The STO phase studied in this work differs from that 
studied by Klein et al. [3§]. The STO phase found 
by Klein et al. arises only in an anti-ferromagnetically 
(J < 0) coupled uncompensated SyF (Mi = MgdiS < 
Msd2S — M2), in which the total magnetization for the 
SyF is opposite to that of mo. However, the STO phase 
found in our study appears in the ferromagnetically cou- 
pled SyF with J > and does not require Mi^M2. Fur- 
thermore, we were not able to reproduced the STO phase 
found by Klein et al. for an uncompensated and antipar- 
allel SyF. We checked the effect of an angular dependence 
of the prefactor Pi that take into account the effects of a 
spin accumulation*^. The boundaries of the white region 
will shift noticeably, but we find no qualitative changes. 
The differences with Ref. 39 might be due to other details 
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in handling spin transport. 

Finally, we note that our approach can be readily ex- 
tended from bi-layer to multilayer systems in which each 
layer is exchange-coupled with its neighbouring layers 
(unpublished). This may provide a novel route to ef- 
fectively synchronize a large network of spin torque os- 
cillators. 

In conclusion, we predict that the ferromagnetically 
coupled SyF can be driven into STO states without the 
need of applying magnetic fields. The resulting STO 



states display large angle precession, therefore generat- 
ing a large power output. In addition to dynamical STO 
states, static canted states are also possible in the same 
structure at slightly different applied current densities. 
Our findings may guide the experimental effort towards 
the field-free STO for real applications. 
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